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We illustrate the energy transfer during planetary flybys as a function of time using 
a number of flight mission examples. The energy transfer process is rather more 
complicated than a monotonic increase (or decrease) of energy with time. It exhibits 
temporary maxima and minima with time which then partially moderate before 
the asymptotic condition is obtained. The energy transfer to angular momentum 
is exhibited by an approximate Jacobi constant for the system. We demonstrate 
this with flybys that have shown unexplained behaviors: i) the possible onset of the 
"Pioneer anomaly" with the gravity assist of Pioneer 11 by Saturn to hyperbolic 
orbit (as well as the Pioneer 10 hyperbolic gravity assist by Jupiter) and ii) the 
Earth flyby anomalies of small increases in energy in the geocentric system (Galileo- 
I, NEAR, and Rosetta, in addition discussing the Cassini and Messenger flybys). 
Perhaps some small, as yet unrecognized effect in the energy-transfer process can 
shed light on these anomalies. 
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1 Planetary fly by s 



The use of planetary flybys for gravity assists of spacecraft became of wide 
interest during the 1960's, when the Jet Propulsion Laboratory (JPL) first 
started thinking about what became the "Grand Tours" of the 1970's and 
1980's (the Voyager missions). The concept was to use flybys of the ma- 
jor planets to both modify the direction of the spacecraft and also to add 
to its heliocentric velocity. At the time many found it surprising that en- 
ergy could be transferred to a spacecraft from the orbital-motion angular- 
momentum of a planet about the Sun, despite the fact it had been known 
since the works of Lagrange, Jacobi, and Tisserand on the three-body prob- 
lem jMoulton 1970~l|Danby 1988] , that the energies of comets could be affected 
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by passing near Jupiter 

Even in the simplest, circular restricted 3-body problem |Danby 1988| , it is 
not that the energy of each object is conserved, only that the total energy 
of the entire system is conserved. Flybys can both give kinetic energy to a 
spacecraft (to boost its orbital velocity) and also can take kinetic energy from 
it (to slow it down for an inner body encounter). 

Hohmann developed a powerful analysis tool for gravity-assist navigation, the 
method of patched conies [Hohmann 19 25 Wicsel 1989]. At JPL clearer under- 
standings of gravity assists were obtained from the works of [Minovich~1 963j 
and [Flandro 1963] . In this technique, a two-body (Kepler) orbit about the 
Sun is first kinematically patched onto a two-body hyperbolic orbit about the 
planet. Then the spacecraft proceeds along the planet-centered hyperbolic or- 
bit during the flyby, and finally a new two-body orbit about the Sun is patched 
onto the post-flyby hyperbolic orbit. 

The question arises in such an approach as to where to do the patch; obviously 
somewhere along the hyperbolic asymptotes, but where? This can be addressed 
by considering how large a region of space is controlled by the flyby planet. 
When the region is reduced to a planet-centered sphere of radius r, and the 
distance between the Sun and planet is given by R, a useful sphere of influence 
is given by the Roche limit, r/R = (Stu/Mq) 1 ^ 3 , where m/M & represents 
the ratio of the planet's mass to the Sun's mass. However, the point where 
the errors in the respective two-body orbits about the Sun and planet are 
equal is determined not by the exponent 1/3, but by the exponent 2/5 as 
derived by Lagrange [Wiesel 19 89J. It actually makes little difference which 
exponent is used to define this matching point, as this whole calculation is 
an approximation to the actual dynamics and, in practice, the gravity assist 
trajectory used for space navigation is computed by numerical integration of 



2 In that same period observations that the solar system could evolve [Hills 1970] 
were similarly met with skepticism. 
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the complete equations of motion. 

However, the patched-conic technique is useful for purposes of searching for 
a fuel-conserving gravity-assist trajectory from a large family of possibilities, 
and especially when multi flybys of more than one planet are used, for ex- 
ample in the Galileo mission [Russell 1992] . In [Flandro 1963] and elsewhere 
[Anderson 1997|TVan Allen 2003] the process has been simply and intuitively 
described and we use it here to gain insights into the nature of the flyby 
anomaly (See Figure [TJ) 




Fig. 1. Geometry of a flyby, modified from [Flandro~1 963j. 

The simple vector velocities in the heliocentric system are added to the orbital 
velocity of the planet (taken to be constant). The initial and final velocities 
in the heliocentric system are Vj and vy. The initial and final velocities in the 
planetary system are and v'*. The velocity of the planet in the solar system 
is v p . The change in kinetic energy per unit mass is 

A£ = (v, • v, - v< • v,)/2. (1) 

A little algebra [Flandro 1963|fVan Allen 2003J gives one 

A/C = v p .(v^-va (2) 

Roughly speaking, in the planetary system which rotates anticlockwise, if a 
satellite in the ecliptic comes from inside the planetary orbit, travels behind the 
planet, and then goes around it counter-clockwise, kinetic energy will be added 
to the spacecraft. Contrarily if a satellite comes from inside the planetary orbit, 
travels in front of the planet, and then goes around it clockwise, kinetic energy 
will be taken away from the spacecraft. 
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Of course one is not violating conservation of energy. The energy (and angular- 
momentum change) is absorbed by the planet that is being flown by. However, 
for such a massive body the relatively infinitesimal change is not noticeable. 
Further, there is (in high approximation) a conserved quantity for the space- 
craft in the barycentric system, Jacobi's integral [Moulton 19701 Danby 1 988 



J = -C/2 = (JC + V) + C = (S ) - uz ■ r x v, (3) 



where {V, C,S} are the potential, rotational-potential, and total energies, re- 
spectively, per unit mass, u is the angular velocity of the planet (system) whose 
vector is aligned with z, the unit vector normal to planet's rotational plane. 
Eq. ([3]) is exactly a constant in the circular restricted 3-body problem, and 
shows how kinetic energy can be exchanged with angular momentum during 
a flyby. 

In this paper we are going to discuss two unusual results that are associated 
with planetary flybys. The first is the possible "onset" of the Pioneer anomaly 
when a major planet gravity assist was used to reach hyperbolic orbit. The 
second result is even more surprising, that during Earth flybys at least three 
craft have exhibited a small velocity increase so that the outbound hyperbolic 
orbits in the Earth system have different energies than the inbound orbits. 

In the following sections we will discuss these anomalies from the energy trans- 
fer viewpoint. For the orbits we discuss we obtained information from JPL's 
Solar System Dynamics web siteOD Note that, although this information is 
precise enough to demonstrate the general energy-transfer properties, to glean 
out the small anomalies we are ultimately interested in understanding, more 
precise orbit determinations are necessary, particularly with regard to the Pi- 
oneer anomaly. We conclude with a discussion. (In Appendix A we describe 
how the original trajectory data was obtained and in Appendix B we give 
further details on the accuracy of the calculations used in this paper.) 



2 "Onset" of the Pioneer anomaly with a flyby boost to hyperbolic 
orbit? 



Analysis of the radio tracking data from the Pioneer 10/11 spacecraft 
[Anderson 1998llAnderson 2002j . taken from 3 January 1987 to 22 July 1998 
(40 AU to 70.5 AU) for Pioneer 10 and from 5 January 1987 to 1 October 1990 
(22.4 to 31.7 AU) for Pioneer 11, has consistently indicated the presence of an 
unmodeled, small, constant, Doppler blue shift drift of order (5.99 ± 0.01) x 
10~ 9 Hz/s. After accounting for systematics, this drift can be interpreted as a 

3 http://ssd.jpl.nasa.gov/ 
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constant acceleration of ap = (8.74 ± 1.33) x 1CT 8 cm/s 2 directed towards the 
Sun, or perhaps as a time acceleration of a t = (2.92 ± 0.44) x 10~ 18 s/s 2 . This 
effect has come to be known as the Pioneer anomaly. Although it is suspected 
that there is a systematic origin to this anomaly, none has been unambiguously 
demonstrated. 

As to the orbits of the Pioneers (see |Nieto 2005] for more details), Pioneer 10, 
launched on 3 March 1972 ET (2 March local time) was the first craft launched 
into deep space and the first to reach an outer giant planet, Jupiter, on 4 
Dec. 1973. During its Earth- Jupiter cruise Pioneer 10 was still bound to the 
solar system. With Jupiter encounter, Pioneer 10 reached solar system escape 
velocity. It is now headed in the general direction opposite to the relative 
motion of the solar system in the local interstellar dust cloud. (Figure [2] shows 
the Pioneer 10 and 11 interior solar system orbits.) 




Fig. 2. The Pioneer 10 and 11 orbits in the interior of the solar system. 

Pioneer 11, launched on 6 April 1973 (ET), cruised to Jupiter on an approxi- 
mate heliocentric elliptical orbit. On 2 Dec. 1974 Pioneer 11 reached Jupiter, 
where it underwent the Jupiter gravity assist that sent it back inside the solar 
system to catch up with Saturn on the far side. It was then still on an ellip- 
tical orbit, but a more energetic one. Pioneer 11 reached Saturn on 1 Sept. 
1979. After its encounter with Saturn, Pioneer 11 was on an escape hyperbolic 
orbit. The motion of Pioneer 11 is approximately in the direction of the Sun's 
relative motion in the local interstellar dust cloud (towards the heliopause). 
It is roughly anti-parallel to the direction of Pioneer 10. 

Why this is of interest is that, as early as about 1980 orbit determination 
results [Anderson 1992] began to show an unmodeled acceleration (see Fig- 
ure [3]), what later became known as the Pioneer anomaly [Anderso"n~1 998j 
[Anderson 200 2] . Note that in Figure [3] there is an apparent onset of the 
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Fig. 3. A JPL Orbit Determination Program (ODP) plot of the early unmodeled 
accelerations of Pioneer 10 and Pioneer 11, from about 1981 to 1989 and 1977 to 
1989, respectively. (This figure originally appeared in [Anderson 19 92].) 

anomaly for Pioneer 11 as it passed by Saturn and reached hyperbolic or- 
bit. If further study demonstrates that the onset of the Pioneer 11 anomaly 
was indeed close to the Saturn encounter time and initiation of the unbound 
hyperbolic orbit, then the obvious question would be if a similar thing hap- 
pened to Pioneer 10 during its Jupiter encounter, when hyperbolic orbit was 
obtained [Nieto 2005j . 

Such further study would be a part of the overall effort to understand the 
time-dependences and precise direction of the anomaly |Nieto 2005] . Since 
the early trajectory data has been retrieved |Turyshev 2006 ] , such work will 
soon be underway. In the meantime we here investigate the mechanism of 
the energy-transfer process. In preparation, in Table [1] we show the orbital 
elements of the Pioneers at the flybys, which resulted in hyperbolic escape 
trajectories from the solar system. 



2.1 Pioneer 11 at Saturn 



As previously stated, it technically is imprecise to consider the total energy of 
the Pioneer 11 spacecraft separately from that of the rest of the solar system. 
Nonetheless, consider the 4-body problem with i) the solar-system barycenter 
(differing from the center because of Jupiter) as the origin of inertial coordi- 
nates and with ii) the potential energy given by the Sun, Saturn including its 
leading multipoles (up to octapole), and Titan. One can then ask what is the 
total energy of the spacecraft, E P , with time. 

In the Saturn-centered frame one expects everything to be symmetric before 
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Table 1 

The osculating orbital elements for Pioneer 11 at Saturn closest approach and Pi- 
oneer 10 at Jupiter closest approach. The quantities are defined in the text except 
for b, the impact parameter, Q, the angle of deflection, (T,d) the time and day 
of closest approach in UTC, and i, the inclination. For Saturn GMg = 37940586 
km 3 /s 2 and for Jupiter GMj = 126712766 km 3 /s 2 . 



OnR.ntitv 


Pio-11 at Saturn 


Pio-1 at .Timitfif 


[km /si 


7.487 


8.477 


vp [km/s] 


31.808 


36.351 


b [km] 


336,855 


869,533 


A [km] 


19,062 


131,279 


e 


1.1176 


1.1150 


G [degrees] 


126.95 


127.49 


i [degrees] 


7.25 


13.77 


T 


16:28:33.370 


02:25:11.135 


d 


01.09.79 


04.12.73 



and after encounter. Going to the solar system barycenter frame, to first ap- 
proximation one a priori might anticipate a continuous, monotonic transfer 
of energy with time. There would be a steady monotonic increase of energy 
before encounter and then a smooth still-monotonically increasing transition 
to a new constant energy after encounter. 

While at Saturn, the time histories for Pioneer ll's kinetic energy, Saturn's 
contribution to the potential energy, the total energy, and the determined 
value of the Jacobi "constant" (all per unit mass) are given in Figure l4lP1 

Figure H] might surprise even though, as we come to below, it was to be ex- 
pected. Its up-down asymmetry is huge and there is nothing like left-right 
symmetry. Pioneer 11 reached a state of positive total energy about 2-1/2 
hours before closest approach to Saturn. But then the spacecraft first gained 
more than the final energy-shift from the planet and then lost some. Does this 
make sense? If we go back to Section [T] it does. 

Reconsider Figured] and Eqs. ([1]) and (jSJ). These last are asymptotic formulae. 
But we can make them a time- dependent formulae by changing Vj to v'(t). 



4 The potential from the Sun varied much more slowly and basically was a bias of 
-94.5 (km/s) 2 . 
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Fig. 4. Energies per unit mass of the Pioneer 11 spacecraft in (km/s) 2 plotted vs. 
time about closest approach to Saturn in hours: the kinetic energy (/Cp = Kp/mp) 
and Saturn's contribution to the potential energy and next the total energy 
(£p = Ep/nip) and the determined value of the Jacobi "constant" (J) computed 
using Eq. (J3j) are given. The latter shows the small residuals about the fit. 

Then, we have 



AK(t) 



(At) 



v 



(4) 



The second term is a constant. By looking at Figure H] one can see that the 
first term gives the spike after perigee with respect to the planet being flown 
by. As the satellite goes around the back side of the planet, v'(t) starts to 
align with v p . The maximum of A/C is reached when v'(t) is parallel with v p . 
This occurs just after perigee. Then as the satellite swings around further, 
v'(£) goes out of alignment with v p so A/C(£) decreases some. The first three 
graphs of Figure H] now makes sense and it is clear what would happen if the 
orbit went the other way (gravity-decrease of total energy). 

The last graph of Figure H] shows that the energy comes from angular momen- 
tum (£). (See Eq. [3j) £{t) has been calculated as a function of time and is 
shown in the third graph of Figure HI z, r, and v can be calculated. Therefore, 
using J, the period r = 2tt/u!, and the declination and right ascension (S z 
and a z ) of the rotation axis z as parameters, a singular value decomposition 
(SVD) fit was done to compare 

£(t) = J-C(t). (5) 
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The three degree of freedom fit results are: J = —121.58 (km/s) 2 ; r = 23.97 
yr; 5 Z = 54.68°; and a z = 268.75°. The resulting period differs from Saturn's 
orbital period of 29.4475 yr by about 19%, but it represents a best fit to the 
total energy in the flyby. The direction of the angular velocity is close to the 
ecliptic pole at a = 270°, 5 = 66.5607°. 



The small deviation of J from a constant is shown in the last graph of Figure 
HI The fact that J is almost the same before and after flyby constitutes Tis- 
serand's criterion for the identification of comets. What matters is that the 
absolute values of J(t), which represent a best fit of the fitting model of Eq. 
(JSJ) to the total energy £ (t) over the flyby interval shown. The time- averaged 
values of J over the flyby intervals are given in the text. Further, the implied 
plot of £ is indistinguishable from the third graph in Figure H] on that scale. 



J is not exactly constant because the real problem is not derived from a static 
rotating potential of the circular restricted 3-body problem. The presence 
of Titan, as well as an oblate Saturn whose rotational axis is not aligned 
with its inertial orbit, both introduce time explicitly into the rotating system. 
This is demonstrated by the deviations from constancy being closest near 
periapsis with respect to Saturn. Note that the Saturn rotational pole is closer 
in direction to the Earth's pole than it is to the pole of the ecliptic, or the 
pole of Saturn's orbit. 



Here and later note that, since the flyby planet is not strictly in a circular 
orbit, the circular restricted three-body problem is only an approximation 
even with no other bodies in the system. The approximation could be made 
better by introducing the elliptical restricted three-body problem. However, 
even with this single improvement the Jacobi integral is no longer a constant 
of the motion. Instead it has a mean value plus time-varying periodic terms 
[Delva 19 79], The constant angular rate u> in Eq. is replaced by a time- 
varying angular rate given by u>(t) = l{t)/r{t) 2 . (Here l(t) is the orbital angular 
momentum per unit mass for the planet about the solar-system barycenter 
and r(t) is the planet's distance from the barycenter.) To the first order in the 
eccentricity e, u(t) = Uq(1 + 2ecosu;o^), with the time t measured from the 
flyby planet's perihelion. However, time variations in the Jacobi integral are of 
concern here only over the time interval of the flyby, which is a small fraction 
of the total orbital period of the flyby planet. The time variation caused by 
the planet's orbital eccentricity is a minor contribution when compared to 
the variations caused by other bodies in the system and by the non-spherical 
nature of the external gravitational potential of the flyby planet (see Appendix 
B). 
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2.2 Pioneer 10 at Jupiter 



Similarly, in Figure [5] we describe the energy-transfer process of the Pioneer 
10 flyby of Jupiter, with the solar-system barycenter as the origin of inertial 
coordinates. The potential energy is given by the Sun, Jupiter with its leading 
multipoles (up to octapole), and the four Galilean satellites. 

Fig. [5] first shows, as a function of time, the kinetic and the Jupiter-caused 
potential energies (per unit mass) of Pioneer 10 during its Jupiter flyby0 The 
third graph in Fig. shows the total energy as a function of time. It is similar 
to the graph in Figure H] except that the total energy exhibits a slight dip 
before its large rise as it approaches periapsis. 
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Fig. 5. Energies per unit mass of the Pioneer 10 spacecraft in (km/s) 2 plotted vs. 
time about closest approach to Jupiter in hours: the kinetic energy (/Cp = Kp/mp) 
and Jupiter's contribution to the potential energy and next the total energy 
(£p = Ep/mp) and the determined value of the Jacobi "constant" (J) computed 
using Eq. ([3]) are given. The latter shows the small residuals about the fit. 

Again using J, r, S z , and a z as parameters, an SVD fit using three degrees of 
freedom was done to compare the two sides of Eq. [5j The results are reasonable: 



J = -234.93 (km/s 



,2. 



10.70 yr; 5 Z = 72.57°; and a z = 297.2°. The 



resulting period is within 10% of Jupiter's orbital period of 11.8626 yr as is 
the direction of the angular velocity with respect to the ecliptic pole direction 



5 The potential from the Sun varied much more slowly and basically was a bias of 
-175.5 (km/s) 2 . 
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at a = 270°, 5 = 66.5607°. 



One therefore again sees that the energy comes from angular momentum. 
Using the values obtained for £{t) and C(t) yields a u J(t), v given by the last 
graph in Figure O Again the Jacobi constant J is still relatively "constant." It 
is not exactly a constant because the real problem is not derived from a static 
rotating potential. The presence of the four Galilean satellites and an oblate 
Jupiter not aligned with its inertial orbit introduces the time explicitly in the 
rotating system. Note that the Jupiter pole is aligned closely to the pole of 
the ecliptic and to the pole of Jupiter's orbit. 



3 The anomalous Earth fly by s 

Earth flybys can be an effective technique for increasing or decreasing a space- 
craft's heliocentric orbital velocity far beyond the capability of its propulsion 
system |D'Amario 1992(Wiesel 1989] . It can also be done repeatedly on the 
time scale of a year. This technique, called Earth Gravity Assist (EGA), de- 
pends on the bending of the geocentric trajectory during the flyby, which in 
turn results in a change in the direction of the spacecraft's geocentric veloc- 
ity vector. This change in direction can cause either a decrease or increase in 
heliocentric orbital energy, depending on whether the spacecraft encounters 
Earth on the leading or trailing side of its orbital path. 

During the flyby the total energy and angular momentum of the solar system 
are conserved. Further, independent of the heliocentric energy change of the 
craft itself, the spacecraft's total geocentric orbital energy per unit mass should 
be the same before and after the flyby. The data indicates this is not always 
true. 

Of course, the conservation of the total geocentric energy before and after 
the flyby is exact in the isolated point mass problem. However, the total 
geocentric energy on entering the Earth's sphere of influence as defined in 
Section. [TJ and on leaving it, should be the same, except for terms in the 
geocentric gravitational potential that depend explicitly on the time. Instead, 
for Earth flybys by the Galileo, NEAR, and Rosetta spacecraft, the geocentric 
orbital energies after the closest approach to Earth were noticeably greater 
than the orbital energies before closest approach. Further, the changes were 
much too large for any conventional time-explicit cause, specifically from the 
Earth's longitudinal harmonics or the motions of the Moon and Sun. So far, no 
mechanism, either external or internal to the spacecraft, that could produce 
these observed net changes in orbital energy has been identified. 

Table [2] shows a summary of the orbital information for these three flybys and 
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also for the interesting Cassini and Messenger flybys. The values in the Table 
come mainly from primary sources, and agree to three+ significant figures 
using orbital information from Horizons. (See footnote [3] and Appendix B.) 
Observe that our altitude, A, is the altitude at closest approach above a ref- 
erence sphere with radius given in Appendix B. These are the radii associated 
with the gravity fields and are not necessarily the measured mean equatorial 
radii. The Messenger numbers, and a few others, come primarily from our 
Horizons calculations. 



Table 2 

Orbital and anomalous dynamical parameters of five Earth flybys, specifically the 
osculating elements at closest approach. The quantities are defined in the text except 
for b, the impact parameter, G, the angle of deflection, (T,d) the time and day of 
closest approach in UTC, (a, 5) the right ascension and declination of the Moon 
in degrees at closest approach, and i, the inclination. GM§ is 398600.4415 km 3 /s 2 . 
The Earth radius, Rq, and mass, Mq, are taken as 6,378.1363 km and 5.9636 x 10 24 
kg. 



Quantity 


Galileo (GEGAl) 


NEAR 


Cassini 


Rosetta 


Messenger 


Voo [km/s] 


8.949 


6.851 


16.01 


3.863 


4.056 


vf [km/s] 


13.738 


12.739 


19.03 


10.517 


10.389 


b [km] 


11,261 


12,850 


8,973 


22,680.49 


22,319 


A [km] 


956.053 


532.485 


1,171.505 


1954.303 


2,336.059 


e 


2.4736 


1.8137 


5.8564 


1.312005 


1.13597 


[degrees] 


47.67 


66.92 


19.66 


99.396 


94.7 


i [degrees] 


142.9 


108.0 


25.4 


144.9 


133.1 


T 


20:34:34 


07:22:56 


03:28:26 


22:09:14.12 


19:13:08 


d 


08.12.90 


23.01.98 


18.08.99 


04.03.05 


02.08.05 


m [kg] 


2,497.1 


730.40 


4,612.1 


2,895.2 


1,085.6 


a [degrees] 


163.7 


240.0 


223.7 


269.894 


107.0 


5 [degrees] 


2.975 


-15.37 


-11.16 


-28.185 


27.59 


Avoo [mm/s] 


3.92±0.08 


13.46±0.13 




1.82±0.05 




Avf [mm/s] 


2.56±0.05 


7.24±0.07 


^ |-0.2| 


0.67± 0.02 


O(0) 


A£ [J/kg] 


35.1±0.7 


92.2±0.9 




7.03±0.19 




AE [J] 


87,600±1,700 


67,350±650 




20,350±560 
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3.1 The first Galileo flyby, GEGA1 



The Galileo spacecraft, launched on 18 October 1989, made use of two EGAs 
to propel it on a trajectory to Jupiter. The first, GEGAl, occurred on 8 
December 1990 at a flyby velocity of Vp = 13.7 km/s and at an altitude 
A = 956 km (flyby radius rp = 7, 334 km). 

Soon thereafter, an analysis of radio ranging and Doppler tracking data re- 
vealed that the pre- and post-encounter data could not be reconciled without 
the introduction of a small but significant Avf ~ 2.5 mm/s velocity increase 
at or near the time of closest approach to Earth |Campbell 1991| . Unfortu- 
nately, during the closest approach the spacecraft was moving too fast and at 
too low an altitude for any station of the Deep Space Network (DSN) to track 
it. Hence the time history of the velocity increase was not obtainable. Even 
so, a later detailed analysis [Antreasian 1998] found that, when expressed in 
terms of the orbital hyperbolic excess velocity Voo = 8.949 km/s, the increase 
was Avoo = (3.921 ± 0.078) mm/s. The equivalent orbital energy increase per 
unit mass at flyby was 

AS = AE/m G « (w 00 )(At; 00 ) « (v F )(Av F ) = (35.09 ± 0.70) J/kg. (6) 
(See column 2 of Table [21) 

Now consider GEGAl from the energy transfer point of view, including the 
potential energies from the Earth, Sun, the Moon, Jupiter, and the Earth's 
quadrapole moment. The last three effects, although interesting, are small. The 
Sun's effect is significantly greater, of order 900 (km/s) 2 , but the change during 
flyby was small. (See Section H] for comments on the Sun's and Moon's po- 
tentials.) That left the Earth's potential energy as the main potential change, 
this not being surprising of course. 

Figure [6] first shows, respectively, Galileo's kinetic and potential energy from 
Earth (per unit mass) during flyby. Next in Figure [6] is shown (per unit mass) 
Galileo's total energy and the variation about the fit to the Jacobi constant, 
-1414.6069 (km/s) 2 , during flyby. For this and the other Earth flybys the zero 
horizontal reference lines represent the average quoted values of J. 

The interesting changes are seen in Galileo's kinetic energy and total energy. 
As with Pioneer 11, there is a peak in kinetic energy transfer after periapsis, 
which peak then comes down to the asymptotic value. This peak is transfered 
to the total energy. The peak is smaller in total relative size than Pioneer ll's 
because of the different scales of the two flybys. The fit to J = —1414.6069 
(km/s) 2 is better than the accuracy of the calculations. 
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Fig. 6. In the solar barycenter system during Earth flyby, Galileo's kinetic energy 
and the potential energy from the Earth and next Galileo's total energy and the 
variation of the Jacobi "constant" (J) about -1414.6069 are given, all per unit mass 
and all plotted in units of (km/s) 2 vs. time from periapsis in hours. 

With only GEGA1 available to suggest an anomalous Earth flyby, and with 
no data to characterize the time history of the anomaly, this result was not 
widely reported. (A second flyby, GEGA2, occurred on 8 December 1992, 
but any potential anomaly was masked by the effects of low altitude drag 
[Edwards 1994] Fn 



3.2 The NEAR flyby (NEGA ) 



On 23 January 1998 the NEAR spacecraft flew by Earth (NEGA) at a velocity 
of vp = 12.7 km/ s and at an altitude in the geocentric system A = 532 km. An 
analysis of the tracking data revealed that an orbital energy increase occurred 
in the vicinity of closest approach. This was true even though this flyby gave 
a negative gravity assist, to reach Eros after the farther-out orbital encounter 

6 GEGA2 had an even lower altitude, 303 km, than GEGAl and had a flyby 
velocity of vp = 14.0 km/s. This resulted in a gap in DSN coverage of nearly two 
hours. In anticipation that an anomalous velocity increase might occur, a Tracking 
and Data Relay Satellite (TDRS) was used to track the Galileo spacecraft during 
the DSN blackout [Edwards 1994] . This time the data revealed a decrease of ~ 5.9 
mm/s in velocity, which could be accounted for reasonably well by the introduction 
of an atmospheric drag model and by a single mass point on the Earth's surface 
underneath the spacecraft trajectory. It was concluded [Antreasian 1998j that if 
an anomalous velocity increase comparable to GEGAl were present, it would be 
masked by uncertainties in the fitting model. 
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with Mathilde. 



Here the encounter came from outside the Earth's orbit, and the peak en- 
ergy transfer occurred just before periapsis. Further, the peak was a positive 
transfer even though the final transfer was negative. This emphasizes the im- 
portance of the vector orientation of the two orbits. 

When expressed in terms of the orbital hyperbolic excess velocity Voo = 6.851 
km/s, the increase for NEAR was [Antreasian 1998] Avoo = (13.46 ± 0.13) 
mm/s. The equivalent orbital energy increase was (92.21 ± 0.89) J/kg. (See 
column 3 of Table [21) Figure [7] shows the kinetic energy and the Earth's 
potential energy on NEAR and NEAR's total energy and the variation about 
the fit to the Jacobi constant of J = —1325.5466 (km/s) 2 , all per unit mass 
as a function of time about periapsis. 




Fig. 7. In the Sun's barycenter system, NEAR's kinetic energy and the potential 
energy from the Earth per unit mass during flyby and next NEAR's total energy 
per unit mass during fiyby and the fit about the determined Jacobi "constant" 
(J = —1325.5466) per unit mass during fiyby are given, all plotted in units of 
(km/s) 2 vs. time from periapsis in hours. 

3. 3 Cassini 

The NEAR result increased interest in the then upcoming Earth flyby by the 
Cassini spacecraft on 18 August 1999. The Cassini flyby altitude was 1171 km 
and its velocity of vp = 19.03 km/s were greater than those for either GEGAl 
or NEGA. It is the only Earth flyby considered here that was prograde relative 
to Earth's rotation. (From a point on the Earth it travels overhead from west 
to east and hence its inclination is less than 90°.) 
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Figure [8] shows the kinetic energy and the Earth's potential energy on Cassini 
and Cassini's total energy and the fit about the determined Jacobi constant 
of J = —1175.3385 (km/s) 2 , all per unit mass as a function of time about 
periapsis. The energy transfer is similar to the Pioneer 11/GEGA cases, but 
the peaks are much smaller. 




Fig. 8. In the Sun's barycenter system, Cassini's kinetic energy and the potential 
energy from the Earth per unit mass during fiyby and next Cassini's total energy 
per unit mass and the fit about the determined Jacobi "constant" (J = —1175.3385) 
per unit mass during flyby are given, all plotted in units of (km/s) 2 vs. time from 
periapsis in hours. 



What, then, about any anomalous energy shift? This time any anomalous 
orbital energy increase was smaller, not surprising given the larger distance 
from the Earth [Guman 2002j . Unfortunately, at the epoch of Cassini periapse, 
an explicit velocity increment of order 2.1 mm/s occurred as the result of 
a series of attitude control jet firings [Roth 2006] . effectively masking any 
observation of anomalous energy shift. Further analysis of the attitude jet 
firings will be needed, but for now we take as an upper limit for any anomalous 
energy shift 10% of the velocity increment currently attributed to the attitude 
jet firings , i.e., -0.2 mm/s. Any significance of the negative sign relative to 
the prograde motion awaits further explanation^ 



Soon after, the Stardust spacecraft had an EGA in Jan., 2001. But there were 
thruster problems which masked any signal. 
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3.4 Rosetta 



The next Earth flyby was by Rosetta, on 4 March 2005. It too gave a pos- 
itive signal of an anomalous energy gain. Rosetta was an ESA craft tracked 
primarily through ESA's European Space Operations Centre (ESOC), and 
in part by NASA's DSN. This provided somewhat independent data anal- 
ysis |Morley 2006| for any obtained Av. The results, primarily from Ref. 
|Morley 2006| , are shown in column 5 of Table [2J There was an anomalous 
velocity increase of Av^ = (1.82 ± 0.05) mm/s. 

Figure [9] shows the kinetic energy and the Earth's potential energy on Rosetta 
and Rosetta's total energy and the fit to a Jacobi constant of J = —1274.7442 
(km/s) 2 , all per unit mass as a function of time about periapsis. To calculate 
this figure we used an Horizons orbit (see footnote [3]) with closest approach 
at 22:09 UTC on 4 March 2005 and at an altitude of 1954 km. The energy 
transfer curves are similar in character to those for Pioneer 11 and GEGA1. 




-4-2 2 4 



Fig. 9. In the Sun's barycenter system, Rosetta's kinetic energy and the potential 
energy from the Earth per unit mass during flyby and next Rosetta's total energy 
per unit mass and the fit about the determined value of the Jacobi "constant" 
(J = —1274.7442) per unit mass during flyby are given, all plotted in units of 
(km/s) 2 vs. time from periapsis in hours. 



17 



3.5 Messenger 



Early analysis results are also in hand from the Messenger craft, which had 
an EGA on 7 Aug. 200fjZ] at an altitude of 2526 km and a flyby velocity of 
10.4 km/s. (See column 6 of Table [2j) This, like NEAR, provided a negative 
gravity assist, taking energy away to eventually achieve Mercury's orbit. But 
the energy transfer curves are very different than NEAR's. 

In Figure [TU] we plot energy functions for the Messenger flyby. (Similar to our 
other calculations, the total potential energy includes contributions from the 
Sun, the Earth and its quadrapole moment, the Moon, and Jupiter.) Unlike 
NEAR, where there was first a positive energy transfer spike, here the figures 
are more an up-down reflection of the Pioneer 11 and GEGA1 figures. There 
is a continuous decrease of energy with time to a negative spike, and then a 
relaxation to the final transfer energy. 
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Fig. 10. In the Sun's barycenter system, Messenger's kinetic energy and the potential 
energy from the Earth per unit mass during flyby and next Messenger's total energy 
per unit mass and the fit about the determined value of the Jacobi "constant" 
(J = —1298.5380) per unit mass during flyby are given, all plotted in units of 
(km/s) 2 vs. time from periapsis in hours. 

Note that now the kinetic and total energies decrease until just after periapsis 
and then increase somewhat to their final values. So this time the energy 
transfer is more a mirror of the GEGAl case rather than mimicking the NEAR 
negative gravity assist. 

The preliminary analysis of the flyby has not, to date, shown any significant 



http : //messenger . jhuapl . edu/missiondesignLive/Tables 



evidence for an anomalous energy. 



4 Discussion 



Having described the dynamical situations of these anomalous fly by s, the ques- 
tion arises if there is any insight that might be obtained. We note again that 
the best-fit values of the rotational frequencies obtained for the Jacobi con- 
stants differ bt 10 to 20 % from the physical values. A priori this could in 
principle be due to missing mass, a misorientation of the invariable plane, an 
added force, or a light speed anomaly, all of which could be mapped into one 
another without better knowledge. 

In this light, it is useful to recall another result from the study of the Pioneer 
anomaly, that there are apparent annual and diurnal terms on top of the 
constant anomaly [Anderson 20 02], These are most apparent in clear data 
sets that came when the Pioneers were at large distances from the Sun. In 
Figure [11] we show examples of this type of signal. 



CHASWP Doppler Residuals For Interval III 




Time in Days from 1 2 July 1 992 00:00 UTC 



Fig. 11. First we give ODP 1-day batch-sequential acceleration residuals using the 
Pioneer 10 data set from 1987.0 to 1998.5 [An derson 2002] . Maneuver times are 
indicated by the vertical dashed lines. The annual term is apparent as an overlay to 
the constant baseline. Next we give CHASMP program acceleration residuals from 
23 November 1996 to 23 December 1996 [Anderson 2002] . A clear modeling error is 
represented by the solid diurnal curve. (An annual term maximum is also seen as a 
background.) 



In discussions of the Pioneer anomaly, these annual/diurnal terms should be 
kept in mind. 

Going on to the Earth flybys, since the mass distribution of the Earth and 
the dynamics of EGAs are supposedly well known, people have started to 
ask if there is at least some phenomenological pattern to the EGA anomalies 
[Anderson 200llLammerzahl 20 06] . 
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Note that the simplest Newtonian version of the problem is the hyperbolic 
orbit in the geocentric system: 



a(e 2 - 1) r F (e+ 1) 
1 + e cos 9 1 + e cos 9 ' 

2M ffi Gl 1/2 



rM m Gi 1/2 



v 2 oo + 



r F 



a 



tan(0/2)' 



(7) 
(8) 



In principle three of these parameters, say Voo, e, and r>, determine the tra- 
jectory. Therefore, one might consider looking for phenomenological patterns 
using such variables. 



Another consideration is the largest non-Earth perturbation, which comes 
from the Moon. (Jupiter has a larger potential contribution but it is the po- 
tential gradient (force) that is important.) In Figure [121 we show the potential 
energies from the Moon during the GEGAl, NEAR, Cassini, Rosetta, and 
Messenger flybys. The Moon's positional dependence could be important, as 
it is part of a separate three-body problem (Sun, Earth, Moon), making the 
actual spacecraft motion closer to a restricted four-body problem. Also, in 
Figure [13] we show the potential energies from the Sun during the respective 
Earth flybys. These are large but do not change much during the flybys. 




Fig. 12. In the solar barycenter system, the potential energies (all per unit mass) 
from the Moon during the (in order) GEGAl, NEAR, Cassini, Rosetta, and Mes- 
senger flybys in units of (km/s) 2 vs. time from periapsis in hours. 



The data points currently available on possible flyby anomalies are limited to 
the two Pioneer flybys of a giant planet and to the Earth flybys of Table |5J 
Unfortunately, the flybys of large satellites orbiting Jupiter and Saturn are 
not reliable because of uncertainties in their gravitational fields, as evidenced 
by Ganymede [anderson 2004|Palguta 2006] . In addition, the detection of any 
small flyby anomaly, with respect to the absolute scale of the velocity in the 
flyby trajectory, requires a spacecraft trajectory relatively free of systematic 
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Fig. 13. In the solar barycenter system, the potential energies (all per unit mass)from 
the Sun during the (in order) GEGA1, NEAR, Cassini, Rosetta, and Messenger 
flybys in units of (km/s) 2 vs. time from periapsis in hours. 

accelerations. These accelerations can be generated by onboard systems, by 
non- conservative forces such as atmospheric drag, or by conservative forces 
from uncertainties in the planet's gravitational field. 

Further, the radio Doppler data must be of high quality and referenced to 
ground-based atomic frequency standards by means of coherent two-way track- 
ing. This limits the number of candidates for further flyby tests, whether of the 
Earth or other planets. However, given the importance of the gravity-assist 
technique for the conservation of rocket fuel, more flybys will occur in the 
future. Now that possible anomalies have been detected, it is just a matter of 
waiting for opportunities and then doing the necessary tweaking of the flyby 
orbit, the onboard systems, and the Doppler tracking in order to optimize the 
possibility for another measurement of an effect. 

In summary, we have presented a description and brief discussion of the physics 
of the energy-transfer process that occurs during planetary flybys. We have 
also presented a series of intriguing real-world results associated with flybys 
that belie our current understanding of the underlying physics. It is hoped 
that further study, which we encourage, can reconcile this situation. 
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A Navigation data 

The navigation results utilized in this paper for Pioneer 10 at Jupiter, Pioneer 
11 at Saturn, [Anderson 2002] . and for the Galileo, NEAR, and Cassini Earth 
flybys [Antreasian 1 998J were produced at JPL/Caltech in Pasadena, Califor- 
nia with the Orbit Determination Program (ODP). The results for Messenger 
were produced at KinetX Corporation in Simi Valley California with their 
Mirage software [Williams 2005] . a navigational system developed at JPL and 
therefore profoundly similar to the ODP. 

The results for Rosetta were produced at the European Space Agencys Euro- 
pean Space Operation Center (ESA-ESOC) in Darmstadt Germany. The soft- 
ware is independent of JPLs ODP, although like ODP it is based on a batch 
least squares procedure. It is intended primarily to fit Doppler and ranging 
data from ESAs 35m New Norcia (NNO) station near Perth Australia, a sta- 
tion of ESAs Intermediate Frequency Modulation System (IFMS), but it can 
fit data from DSN tracking stations as well. The Rosetta Radio Science Team, 
which depends on both IFMS and DSN tracking data, is a collaboration of 
ESA and NASA scientists [Patzold 200TallPatzold 2001b] . 

All these software systems depend on similar fitting models |Moyer 2 000 1. 
The dynamical models are referred to the J2000 inertial reference frame of 
the JPL export ephemerides, with Barycentric Dynamical Time (TDB) as the 
independent variable for the spacecrafts vector position and velocity. In ad- 
dition to the Newtonian attraction of the planets and the Moon, including 
a high degree and order truncated Legendre expansion of the Earths gravi- 
tational potential, the models include post-Newtonian corrections consistent 
with the theory of general relativity, models for the reaction forces from solar- 
radiation pressure on the spacecraft, models for orbital trim maneuvers, and 
stochastic small forces from gas leaks, nonisotropic thermal emission, and un- 
coupled control jets. The numerical integration of the equations of motion 
is referred to the Earth's center of mass for the Earth flybys, although the 
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integration is dynamically consistent with the inertial solar-system barycen- 
ter as the dynamical center. The Doppler and ranging data are corrected for 
transponder delay, refraction in the Earths troposphere and ionosphere, and 
for refraction by interplanetary plasma, the latter being of small concern for 
the Earth fly by s. Corrections are also applied for variability in the length of 
day (the Earth's rotation) and for polar wandering in the Earth's body-fixed 
coordinate system. 

The data delivered for analysis depends on the ground system, whether DSN 
or IFMS, and whether closed loop or open loop. For all the Earth flybys, the 
data is available in the closed loop mode and are extracted at the stations by 
counting cycles of a sinusoidal signal recorded by digital receivers as a function 
of time (UTC). The data acquisition and data processing is similar for the DSN 
and IFMS and the delivered data are compatible. The raw data (cycle count) 
for the Earth flybys is available in archival tracking data format (ATDF). The 
first level of data processing is done by the DSN or by IFMS and is delivered as 
an Orbit Data File (ODF). The ODF consists of samples of Doppler frequency 
shift, defined as the difference of cycle count at a predetermined Doppler 
integration time TC divided by TC and referenced to a time variable uplink 
frequency as recorded by the transmitting station |Moyer 2000 . 

The result is called closed-loop two-way Doppler data when the transmitter 
and receiver are at the same station, and three-way Doppler data when the 
transmitter and receiver are at different stations. Ranging data is delivered on 
the ODF at the sample times recorded by the ranging receivers at the station 
and is in ranging units RU, which can be converted to a round-trip light time 
in UTC seconds |Moyer 2000] . Any further data processing is done by the 
data analyst, and it is done with a software system such as JPLs ODP that 
accepts an ODF file as input data. In that sense, the data displayed in Table 
|2l including the anomalous Earth-flyby results, can be taken as processed 
tracking data. There is no need to go back to the processed data on the ODF, 
or the raw data on the ATDF, except for purposes of running independent 
checks on the data processing that leads to the results of Table [2J 



B Flyby calculations 

B. 1 Major planet flybys by the Pioneers 

The Pioneer calculations used the potentials from the Sun, Jupiter, and Saturn 
(GMs), the Earth's effect being very small and hence possible to ignore. For 
the Pioneer 11 flyby of Saturn, the specific Newtonian approximation used to 
calculate £ was, 
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where the velocity v is with respect to the solar-system barycenter (SSB) at 
distance r with solar constant GM Q . 



The gravitational constant for the Saturn system is GM$- Saturn's oblateness 
coefficients are given by , Jf , and Jg for a reference radius Rs, P n is the 
Legendre polynomial of degree n, and the sine of the latitude is obtained by 
the scalar product of the unit vector k directed to Saturn's north pole and the 
unit vector t$ directed from Saturn's barycenter to the spacecraft. The right 
ascension and declination of the Saturn pole in J2000 coordinates is given by 
ctp and Sp so that k = [cosa^ cos 5^, sin cos^p, sin<5p]. The sixth harmonic 
coefficient J$ is negligible for three-place accuracy in the energy calculation, 
but it is included. The values of the Saturn constants used in the calculation 
are given in subsection IB. 31 



Jupiter simply adds a constant bias to the energy over the time interval of 
the flyby, and can be ignored. The flyby is referenced to the barycenter of 
the Saturn system, rather than the center of the planet, with rs the distance 
between the Saturn barycenter and the spacecraft. By using the barycenter, 
the satellites of the system are accounted for, and the approximation to the 
energy is more accurate, just so long as the spacecraft does not approach a 
satellite at close range. 



The energy calculation for the Pioneer 10 flyby of Jupiter is formulated exactly 
the same as in Eq. (IB.lj) for Saturn, but with the subscript /superscript S 
replaced by J. The values of the Jupiter constants are also given in subsection 



According to the Horizons trajectory, the closest approach of Pioneer 11 to 
the Saturn barycenter was 79446 km on 1 September 1979 16:29:23.553 ET or 
16:28:33.370 UTC. The closest approach of Pioneer 10 to the Jupiter barycen- 
ter was 202756 km on 4 December 1973 02:25:55.318 ET or 02:25:11.135 UTC, 
with times referred to the spacecraft. The Pioneer trajectories archived on 
Horizons possibly could be improved for these two flybys by new fits to the 
recently-retrieved early Doppler data |Nieto 20051 Turyshev 2 006 1. 
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B.2 Earth flybys 



The energy calculations for the Earth flybys are also based on spacecraft tra- 
jectories archived in JPL's Horizons system (see footnote As such, the 
accuracy of the spacecraft orbits is limited by whatever trajectory file was de- 
livered to Horizons by a navigation team working for a particular flight project 
at the time of the flyby Presumably, all the delivered trajectories provide a 
good fit to the DSN (or ESOC) radio Doppler and range data generated dur- 
ing the flyby. However, even though we can assume the trajectories provide 
a good fit to the data, the use of the Horizons system to compute orbital 
energies has its limitations. 

A conservative estimate of the error introduced by the Horizons system can 
be obtained by doing the energy calculations both by the "observer" method 
and also by the "vector" method and comparing the results. We conclude that 
the Horizons system introduces an error of no more than 0.1 (km/s) 2 . We 
proceed with the energy calculations with the Horizons' data, but restrict the 
application of £ to three significant digits. 

In the end, any inconsistencies between the inbound and outbound Earth 
flyby data have been reconciled by the introduction of a fictitious maneuver 
at perigee in the direction of the spacecraft motion. For the Galileo, NEAR, 
and Rosetta flybys, such a maneuver definitely is needed in order to fit all the 
inbound and outbound data to the noise level with a single trajectory. 

Another potential problem with the energy calculations is that the total or- 
bital energy per unit mass, £ , is computed with the Newtonian approxima- 
tion in inertial coordinates. However, the spacecraft trajectory and the JPL 
ephemerides are computed consistently to post-Newtonian order (order v 2 /c 2 ). 
It is not solely that post-Newtonian terms of order 10~ 8 should be added to 
the energy calculation, although this could be done, but rather that the con- 
stants and coordinate positions of the planets and the Moon would be different 
if the solar-system data were fit with a Newtonian model as opposed to the 
relativistic model actually used. 

However, if the coordinates are evaluated at the dynamical time TDB of the 
ephemeris, and the constants of the ephemeris are used in the calculation, 
the value of £ should be good dynamically to at least seven significant digits. 
The Horizons trajectories limit the £ calculation, not relativity considerations. 
Consequently, we include only four principal bodies in the Earth flyby calcula- 
tions; the Earth, along with its oblateness coefficient, the Sun, the Moon, and 
Jupiter. This assures that the result for £ is good to three significant digits. 

The Newtonian approximation used to calculate £ for Earth flybys is the 
following: 
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GM M GM 



GM, 



a 



1_ J°(**) P 2 (sin5) 



(B.2) 



where the velocity v is with respect to the SSB, and the gravitational constants 
are GM Q for the Sun, GMm for the Moon, GMj for the Jupiter system, and 
GM @ for the Earth. The Earth's oblateness coefficient is J® for a reference ra- 
dius R 9 , P2 is the Legendre polynomial of degree two, and we approximate the 
true geocentric latitude by the declination 5 in inertial coordinates. Because 
the J® term adds a maximum magnitude of 0.02 (km/s) 2 to £ for the NEAR 
trajectory, it can be approximated by the declination and still be below our 
level of assumed systematic error. The spacecraft-body separation distances 
are r , tm, O, and r ffi for the Sun, Moon, Jupiter's barycenter, and Earth. 

B. 3 Numerical values of constants 

The solar and Saturn constants used in the Pioneer 11 calculation are the 
following (see footnote [3]) : 



GM Q = 132712440018 km 3 /s 2 

GM S = 37940586 km 3 /s 2 = GM /3497.898 

R s = 60330 km 

J| = 0.0162906 

Jf = -0.000936 

JS = 0.000086 

a s = 40.58364 deg 

6° = 83.53804 deg 

The additional Jupiter constants used in the Pioneer 10 calculation are (see 
footnote [3]). 

GMj = 126712766 km 3 /s 2 = GM /1O47.3486 

Rj = 71492 km 

J§ = 0.01469645 

Jf = -0.00058722 

Jjf = 0.00003508 

a J p = 268.0567 deg 

5/ = 64.4953 deg 

The additional constants used in the Earth flyby calculations are (see footnote 



E]): 

GM e = 398600.4415 km 3 /s 2 
GM M = 4902.798 km 3 /s 2 
R m = 6378.1363 km 
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J 2 ® = 0.0010826269 
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